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In cuprates superconductors (as polycrystal and single crystal samples) unusual hysteretic temperature be- 
havior of elastic properties is frequently observed. In this work from the uniform point of view detailed analysis 
this abnormal temperature behavior of elastic properties is given. It is shown that the hysteretic temperature de- 
pendencies of elastic properties of YBa2Cu307_« crystal are strongly anisotropic and are connected to hysteretic 
behavior of the C3333 modulus while the Clin and C2222 moduli have no such dependence. The analysis of the 
elastic constant tensor on the basis of the microscopic model has allowed to draw a conclusion, that the hysteretic 
behavior of the C3333 modulus of YBa2Cu307_4 crystal is caused by temperature dependent renormalization of 
constants of interaction of apex oxygen atoms with copper atoms and are connected apparently to formation of 
charge ordering. 



I. Introduction 

Already the first investigations of high temperature su- 
perconductors 01211110] have revealed a variety of essen- 
tial features in the temperature dependencies of velocity 
(v(T)) and attenuation (a(T)) of ultrasound. So for the 
most well investigated isostructural series of HTSC com- 
pounds (RE)BaaCu307_i5 (RE - rare earth atom) in the 
temperatures region 190 - 230 K step changes of v(T) 
have been observed |S1 [B] , indicating on lattice instabil- 
ities and/or phase transitions at temperatures previous 
to the superconducting transition. Surprising there was 
that at heating and cooling of a sample the tempera- 
ture dependencies of velocity (attenuation) did not co- 
incide. A distinct thermal hysteresis was observed be- 
tween 60 and 230 K E] . Further investigations have 
shown that the hysteretic behavior of w(T)is observed 
for all isostructural series (RE)Ba2Cu307_5 03 ^J, and 
also for the Bi [12 EH El 021 , Pb 0)], Tl Q7] HTSC 
compounds. Also appeared that for polycrystal sam- 
ples the value of the hysteresis depends on structure of 
polycrystal samples. So hysteretic behavior of u(T)was 
observed in textured polycrystal samples [E! and poly- 
crystal samples with relatively large crystal grains (^50 
/zm) 05]. I n nne grinding, dense polycrystal samples [20] 
and also polycrystal samples obtained by MPMG (melt- 
powder-melt-growth) method |21j such features were not 
observed. The results marked above specified a possi- 
ble connection of anomalous temperature dependencies 
v (T) with physical processes proceeding on crystal grain 
boundaries and, hence, for single crystals hysteretic be- 
havior should not be observed. However, further inves- 
tigations have revealed the hysteretic behavior of v(T) 
in (RE)Ba 2 Cu 3 07_ 5 single crystals [23 E3) So in [22] 
for ultrasound with frequency 10 MHz the hysteretic be- 



havior of the C3333 and C1212 elastic constants was ob- 
served , while hysteretic behavior of Cnn, C2222, C1313 
and C1212 was not observed. When studying of elastic 
properties of single crystals YBa 2 Cu307_,$ at 100 Hz 
hysteretic behavior of both elastic and torsion modulus 
were observed [2B]- This fact disagrees with results of 
[22|, where the hysteretic behavior of the C1313 and C1212 
modules was not observed (module of torsion is deter- 
mined by these elastic constants). Observation in single 
crystals of hysteretic behavior of w(T)puts under doubt 
the models offered before for an explanation of these de- 
pendencies (see, for example, EH])- Really, it is gen- 
erally believed that in cuprate superconductors anoma- 
lous temperature behavior of ultrasonic wave velocity is 
result from redistribution of oxygen atoms in the Cu- 
O chains or/and phase transitions in twinning area of 
crystals [21 125) . In this connection it is necessary to 
note, that the hysteretic behavior of velocity of ultra- 
sound also was observed in (RE)Ba2Cu40s polycrystals 
PB1 ffi\ . This material is stable down to 850 °C and 
crystals are untwined. Besides while investigating of a 
linear expansion coefficient of YBa2Cu307_5 single crys- 
tals in 28 was shown that in chains disordering processes 
of oxygen atoms proceed at temperatures above 300 °C. 
Let's pay also attention that the hysteretic behavior of 
v(T) recently was observed by us and in Ba^Ki-jRiC^ 
compound 29 . For Ba^Ki-^BiOs and (RE)Ba 2 Cu40g 
compounds the hysteretic behavior of ultrasound can not 
be closely related to redistribution of oxygen atoms or 
twinning. Earlier we offered a phenomenological model 
describing hysteretic temperature behavior of elastic con- 
stants of HTSC single crystals [23] • However the micro- 
scopic mechanism resulting in so anomalous temperature 
behavior of velocity of ultrasound and essential distinc- 
tion in elastic properties of polycrystal and single crystal 
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samples of HTSC compounds remains obscure till now. 
The solution of this task is important for the following 
reasons. Recently for an explanation of the phase dia- 
gram of cuprate superconductors (remaining mysterious 
after 13 years of research [2]) the local micro-strain of 
Cu-0 bonds and charge ordering (formation of striped 
phases) jHJ were used. The micro-strain of Cu-0 bonds 
should result in appreciable changes of elastic properties 
of crystal. Besides, hysteretic behavior of wave veloc- 
ity of ultrasound was discovered recently in MgB2 com- 
pound |S21- The huge temperature hysteresis of veloc- 
ity of a sound and internal friction was observed also in 
La .8Sr .2MnO 3 single crystals [33]. 

In the given work is shown that the hysteretic behav- 
ior of velocity of ultrasound, as in single crystal, and in 
polycrystal samples of coprate superconductors may be 
explained from the uniform point of view. In a final part 
the possible microscopic mechanism of hysteretic behav- 
ior of ultrasound is discussed. 

II. Anomalous temperature dependencies 
of elastic properties of single crystals 

To develop a detailed microscopic model of hysteretic be- 
havior of ultrasonic velocity and to analyze temperature 
dependencies of ultrasonic wave velocity in polycrystal 
materials we must know the components of elastic tensor 
which have anomalous temperature dependence. How- 
ever to the present time there are only two experimental 
works devoted to detailed analysis of temperature de- 
pendencies of elastic properties of YBa2Cu307_^ single 
crystals but as I have marked in introduction, from [22] 
and [53J directly it is impossible to make an unambiguous 
inferences about the temperature dependencies of elastic 
modules. First of all I show that the experimental results 
of [22] are consistent with results |23] quantitatively in 
the assumption, that anomalous temperature dependence 
the elastic modulus C3333 has only. 

Let's calculate resonant frequencies of a sample at de- 
formation of a torsion and compression in view of inertia 
of cross movement of elements of a sample. We will con- 
sider the case when the laboratory system of coordinates 
is oriented arbitrarily in reference to the crystallographic 
system of coordinates and the calculations are feasible by 
a variation method. 



A. Orientations m_Lc, fc||c; fc_Lc, u\\c 

For the given experimentally investigated in [23] orien- 
tations of laboratory and crystallographic systems of co- 
ordinates (u is the displacement vector of an ultrasonic 
wave; c is the unit vector of the crystallographic direction 
[001]; k is the wave vector) in a sample the deformation of 
a torsion and a bend is created. For evaluation of natural 
frequencies of oscillations of a sample as generalized coor- 
dinates it is convenient to pick the tensor fly = ^ ( ^ i /g x . 



(see, for example, [2H]), which is closely related to the de- 
formation tensor (uy ) f2y = |36| . If the tensor 
Ojj is known, that the deformation may be restored to 
within moving a sample as the whole |37| . When measur- 
ing natural frequencies of a sample value f2y is usually 
small and constant. For small values of Ojj we can ex- 
pand a vector and tensor of deformation into a Taylor 
series of f2y and take into account linear terms only 



Ui = Bijk(xi,X2,xa)Cljk', 
Uij = Aijki{xi,x 2 ,x 3 )flkl, 



(1) 



Then the free energy of the deformed sample may be 
written as follows 
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r = - J c ijkl u ijUkl dV 

V 



Akl ra dV I Q mn £l r q 7 (2) 



where Cy/y is a tensor of elastic modules. When writing 
down kinetic energy as 

T= - J pu 2 dV = i J pBijk&jkBimiQmidV (3) 

V V 

we obtain the following expression for a Lagrangian 
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L T F QPBjkmlQjkQml 2 Csgrp^sq^rp: (4) 



where 



Bjkml = J BijkBimidV; 
V 

Cmnrn — J Ci^kl^-ijrLrnA-klrnd^ • 
V 



(5) 



The Lagrangian (4) leads to the equation of motion for 

pBijral^ljnl ~\- Cij r pQ r p — 0. (6) 

Now using the equations of motion (5) we may ob- 
tain the characteristic combined equations for definition 
of natural frequencies of oscillations of a sample 
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For the analysis of experimental results [23] we shall 
choose a sample in the shape of the elliptic cylinder. It 
will allow to take into account all characteristic dimen- 
sions of experimentally studied samples |23j and to derive 
an expression for f2y in the closed form, that essentially 
simplifies evaluation of resonant frequencies. Generally 
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the stress tensor is determined by the following equations 
and boundary conditions 36 : 



dtJij 
d.r k . 



0. 



(TijUj = p t , (c) 



(8) 



where rij is the normal unit vector to a surface of a sam- 
ple, pj is the force acting on unit area of a sample. For 
the case of an isotropic elliptic cylinder the er^ stress 
tensor is a linear function of coordinates [12] and the 
equations (7b) are satisfied identically. Thus, the stress 
tensor is determined only by the equations (7a), which, 
generally speaking, have the same form as for isotropic, 
and an anisotropic material. Hence, the tensor of strains 
and for the anisotropic elliptic cylinder will be a linear 
function of coordinates and may be found on the basis of 
the analysis of boundary conditions and has the following 
form 
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where M is the torsional moment; Sijki is the tensor 
of elastic compliances in laboratory system of coordi- 
nates and is related to the tensor in crystallographic 
system through the transformation matrix as Si>j>k'i' — 
A^A^A^AySijki By using (1) the tensor fi^ may be 
written as 
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In the case of u_Lc, k\\c the equation of motion for Qij 
(6) leads to the following expression for natural frequen- 
cies of sample oscillations 

1 C1313C2323 Ylab 



oj 2 = 



P L2 [(-a) Cl313 
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b 2 ' 



(11) 



where p is the density of a sample, L is the length of a 
sample along a torsion axis; Ciju is the tensor of elastic 
modules expressed in laboratory system of coordinates. 
In the case fc_Lc, u\ |c the characteristic equation for res- 
onant frequencies of a sample have been obtained in a 
similar way and because of complexity has been solved 
by numerical methods. 



B. Orientations 



u\\c,i 



u\ \c, k-Lc 



For the given orientation of laboratory system of co- 
ordinates longitudinal vibrations are raised in a sample. 



Calculation of natural frequencies of longitudinal vibra- 
tions of a sample was performed on the basis of the equa- 
tion of motion for a vector of strains: 

d 2 u 3 CijkiUijUki ( d 2 u 3 



dt 2 
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(12) 



where the ETy tensor is determined on the basis of the 
condition ity = -IIyit33 ; i?n = Uf 3 - ILui'&n = #21 = 
II13II23 — II12; ??22 = n| 3 — II22; 5a is the cross-sectional 
area of a sample. The equation (11) is obtained by a vari- 
ation method and for isotropic materials coincides with 
the Relay equation The harmonic solution of the 

equation (11) can be deduced by the method of parti- 
tioning variables and results in the following expression 
for natural frequencies of oscillations of a sample 



CijkiRijTlki 



P 



n 2 13 + n 2 3 + ( 2 -^y^ 



(13) 

where n is an integer. 

Let's analyze temperature dependencies of the natu- 
ral frequencies w(T), calculated on the basis of the ex- 
pressions obtained above . First of all, one additional 
remark need to be made. Recently we have shown 30 |, 
that the C3333 (T) temperature dependencies observed in 
|22| may be qualitatively described within the framework 
of model of a strongly correlated bistable sublattice and, 
hence, simulation of natural frequencies oj(T) of a sample 
also may be performed within the framework of the given 
model. However here I shall act as follows. At calculation 
of temperature dependencies of natural frequencies it is 
convenient to choose as a modelling temperature depen- 
dence of 6*3333 (T) the C3333(T) dependence measured in 
[22] > that enable not only to simulate temperature depen- 
dencies but in details allows to compare results of works 
[221 and [231 • The calculated temperature dependencies 
w(T) are depicted in figure 1. For all investigated orien- 
tations calculated and experimental dependencies w(T) 
are in excellent agreement. I pay attention, that for the 
orientation u\\c, k _L c a strong temperature dependence 
of w(T) is caused by excitation of bending oscillations of 
a sample which natural frequency is determined by the 
C3333(T) module . 

III. Analysis of abnormal temperature 
dependencies of ultrasonic velocity in 
polycrystal samples 

For calculation of effective values of a tensor of elastic 
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Figure 1. Temperature dependencies of the relative resonant frequency (w(T) — oj (To)) /to {To) of oscillations of a 
single crystal sample for various orientations of the wave vector k and the vector of displacement u of a ultrasonic 
wave . Lines represent simulations; points are experimental results |23|: the direction of change of temperature is 
shown by arrows 



modules of polycrystal materials, it is necessary to know 
topological type of the heterogeneous medium. From 
topological point of view all heterogeneous mediums may 
be subdivided on matrix and statistical mediums. In ma- 
trix mediums the main phase envelops all other phases 
and forms connected space at any volume concentration 
of the main phase. In statistical medium all phases are 
distributed by a casual fashion not forming regular struc- 
tures. For definition of topological type of the polycrystal 
medium we shall break volume of a sample on elemen- 
tary cubic regions (domains) with the characteristic size 
Lq , equal to the average size of a crystal grain. Let 
the crystal grain has orthorhombic symmetry. Each sep- 
arate crystal grain is oriented arbitrarily relative to the 
chosen axis x 3 determining propagation of an ultrasonic 
wave. All set of domains of dissection breaks up to three 
subsets in which domains with orientation of crystallo- 



graphic axes of a crystal grain (a, b, c) belonging to the 
intervals: 

Intr(a) E jo < L(ax 3 ) < j; ^ < l(ax 3 ) < 2tt| ; 
Intr(c) E jo < l(cx 3 ) <^^- ~ 2lT ] 5 

Intr(b) £ jo < l(bx 3 ) < J; ^ < l(bx 3 ) < 2tt| . (14) 

The given splitting can be performed always as it is 
possible to show that for Euler angles the given relation 
is always fulfilled: 

(f<»<o)*( v >:,*>J). da) 

It is obvious, that for the chosen subset of domains 
the distribution function of orientations of crystal grains 
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f(6) has a maximum in the intervals (13). And, texture 
of a sample is determined by a relation of numbers of fill- 
ing of the subsets defined above. Further, for the given 
subset of domains we replace the tensor of elastic modules 
of each domain with average value of the tensor of clastic 
modules < Ciju > v (r] = Intr(a), Intr(b), Intr(c)) [3~9"| 

1 



{Cijki) 



8tt 2 



C m (9, <p, 1>) f (0, <p, i>) sm(0)dcpdi>de. (16) 



When performing concrete calculations the distribution 
function f(6,if>,(p) was chosen in the form f{d,ip,(p) = 
[cos(6*)] 2 ™ [39], n - parameter (at n — > or n — * oo we 
have random or textured distribution of orientations of 
crystal grains, respectively). As a result, polycrystal ma- 
terial can be considered as the three-phase heterogeneous 
medium. Topological type of medium is determined by a 
relation of the relative volumes Vg , Vg ,Vg and the perco- 
lation threshold X c ( X c = 0.307-^0.325 for the site prob- 
lem on a cubic lattice @HJ). So the conditions for matrix 
and statistical mediums have the form Vg < X c (matrix 
medium concerning phases Vg ,Vg) and Vg , Vg ,Vg > X c , 
respectively. 

C. Finely grained dense polycrystal ceramics 

For YBa2Cu307_a compound at S » crystal grains 
has orthorhombic symmetry. Orientation of crystal 
grains is arbitrary for a nontextured sample and for rel- 
ative volumes of separate phases we have the following 
condition Vg w Vg f» Vg w ^ . The calculation of average 
values of a tensor of elastic modules has shown, that for 
compound YF^CuaOy-a we have the following relation 
< Cnii >~< C2222 >• Hence, in this case polycrystal 
material is matrix medium and the matrix is formed by 
crystal grains with relative volume Vg + Vg > X c . For the 
matrix material the effective elastic modulus may be cal- 
culated using the virial approach based on decomposition 
of an effective tensor of elastic modules in a series on rel- 
ative volume of a phase Vq < X c and in linear approach 
expression for the effective module looks like 



C* = (6 



1111,2222 



1 



3((C 33 33>-(C; 



1111,2222 



)) 



2 (63333) 4- (61111,2222) 



Vg 



(17) 

Using expression (14) and average values of a tensor of 
elastic modules (15) we have calculated temperature de- 
pendence of ultrasonic velocity for dense finely grained 
polycrystal materials. At calculation was considered, 
that the temperature dependencies of elastic modules of 
separate crystal grains and single crystals are identical. 
Results of calculation are depicted in figure 2. The main 
feature of temperature dependencies in this case is ab- 
sence of abnormal behavior of temperature dependencies 



of ultrasonic velocity though the elastic modules C3333 of 
separate crystal grains had pronounced hysteretic tem- 
perature dependence. Let's pay attention, that the calcu- 
lated temperature dependencies are in a good agreements 
with experimental ones |2()| . 

Thus, for not textured finely grained polycrystal ma- 
terials the temperature dependencies of velocity of ultra- 
sound are determined by average values of the < Cxin > 
, < C2222 > crystalline elastic modules. 

D. Textured polycrystal materials 



In the case of textured polycrystal material for the rel- 
ative volume occupied by the crystal grains oriented, for 
example, in c direction the following condition is fulfilled 
Vg > Vg, Vg and the relative volumes Vg, Vg are much less 
than the percolation threshold X c . In this case we have 
the matrix topology of polycrystal medium, but the ma- 
trix is formed by crystal grains with average value of the 
elastic modulus C3333 . For calculation of temperature 
dependencies expression (14) can be used in which the 
following replacement < C3333 ><^< 61111,62222 > have 
been performed. Calculated temperature dependencies of 
ultrasonic velocity for textured polycrystal samples with 
relative volume V c « 0.9 are depicted in figure 2b. In 
this case hysteretic behavior of w(T)is observed. Similar 
temperature dependencies were observed experimentally 
for textured samples in |18j . 

E. Materials with topology of a statistical 
mixture 



Let's consider polycrystal materials, when for relative 
volumes of separate sets of crystal grains the following 
condition is satisfied Vg » 0.5 > X c ; Vg U Vg » 0.5 > X c . 
In general case for calculation of an effective tensor of 
elastic moduli detailed distribution of inhomogeneities in 
a sample must be known 35 . However temperature de- 
pendence of ultrasonic wave velocity can be analyzed on 
the basis of the percolation theory. Really, from the per- 
colation theory follows immediately 001 that the whole 
volume of a sample of a statistical mixture can be divided 
into two connected clusters with different orientation 
of crystal grains Intr(c) and Intr(a), Intr(b), respec- 
tively. And, apparently, that properties of a statistical 
mixture should be symmetrical concerning replacement 
< 63333 >^< 61111,62222 >■ For a wave length of ultra- 
sound the following relation X sound < L duster w L samp i e 
is fulfilled and elastic waves propagate in both clusters. 
When an elastic plane wave propagates in nonuniform 
material wave front will be deformed and the spatial mod- 
ulation of amplitude will appear. It will take place even 
in a case when end faces of a sample are absolutely par- 
allel and flat. The value of 
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Figure 2. Temperature dependencies of the relative velocity of ultrasound ((v(T) — v(To))/v(Tq) for a,b and 
V(T)/v(To) for c) in Y BC12CU3O7-S polycrystal samples. Lines represent simulations; points are experimental 
results: a - for dense, small crystal grains sample |20):b - textured sample |18) : c - sample with relatively large 
crystal grains (~50/xm) |19| : the direction of change of temperature is shown by arrows 



a signal received by a detector sensitive to amplitude of 
strains will be proportional to a wave amplitude averaged 
over the detector surface 



velocity in a case of a statistical mixture can be obtained 



u Q = k {A) s ; (A), 



~ J Y^ A i( x ,y)dxdy, 



(18) 



where summation is carried out over all connected clus- 
ters. Elastic waves propagating in connected clusters are 
detected as a single wave with effective half- width and ef- 
fective standing of a maximum < t > . The average value 
< t > is determined by a relationship of amplitudes Ai 
and times ti of elastic waves propagating in each clus- 
ter and in a linear approach may be represented as the 
ensemble average 



(19) 



where 7 is the parameter of the given detector sys- 
tem. Expressing the amplitude of a passed wave as 
Ai = Aq exp (— (XiV) the following expression for average 



T I ^ 



[1 + exp (aiL 1 - a 2 L 2 ) 



(1$f}) ex p b ( a i L i - "2^2)] 



where Vi , a, are the velocity and the decay factor of an ul- 
trasound wave, Li is the average size of a cluster; So is the 
average sectional area of a cluster; S is the effective area 
of a detector. As follows from equation (19) , < v(T) > 
temperature dependence is being determined by v\ (T) or 
v 2 {T) temperature dependence depending on a relation 
of a.i(T) decay factors. I have simulated < v(T) > tem- 
perature dependencies of ultrasonic wave velocity taking 
into account experimentally measured temperature de- 
pendencies of a decay factor |19II22| (see Figure 2). Tem- 
perature dependence at high temperatures is determined 
by temperature dependence of ultrasonic wave velocity in 
the cluster with orientation of crystal grains Intr(c). At 
decreasing of temperature the decay factor of ultrasound 
in the cluster with orientation of crystal grains Intr(c) is 
sharply increasing while in the 
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Figure 3. Temperature dependence of incoherent lattice fluctuation in YBa2Cu307_a crystals, 
imental results [45] . Lines are calculated temperature dependencies: 1 - Debye model; 2,3 
correlated sublattice without and with taking into account disordering of oxygen atoms in Cu(l)-0(1) chains. On 
the inset - expansivity versus temperature of untwined YBa2Cu307_^ single crystals for S ~ 6.95 (curve 1) and 
S ~ 7(curve 2) |2^]. 



cluster with orientation Intr(a), Intr(b) is decreasing 
[TT?I |221 an d as a result velocity of ultrasound will be 
determined by the cluster of the type Intr(a), Intr(b). 
The given temperature dependence of u(T) represents in 
some sense a combination of the temperature dependen- 
cies considered above for textured and finely grained ma- 
terials. It is interesting to note, that a similar sort of 
temperature dependencies were observed experimentally 
in samples with large crystal grains (50 micron) (see 
fig. 2). Thus, quantitative analysis of experimentally 
observed hysteretic temperature dependencies of ultra- 
sonic wave velocity in YE^CuaC^-a polycrystal mate- 
rials (finely grained, textured, macrocrystalline) at fre- 
quencies from 100 kHz up to 20 MGz has shown, that 
as well as in case of single crystal samples, the hysteretic 
behavior is determined by hysteretic temperature behav- 
ior of the modulus C3333 while the other elements of a 
tensor of elastic moduli have no abnormal temperature 
dependence. 



F. On microscopic mechanism of hysteretic 
behavior 

The fact of existence of strong anisotropy of hysteretic 
behavior of elastic modulus is important from the point 
of view of clearing up of the microscopic mechanism of 
so unusual behavior. Really, the essential contribution 



to anisotropy of a tensor of clastic modulus is given by 
interactions of atoms of the nearest coordination spheres 
of crystal. Having expressed a tensor of elastic moduli 
as a function of constants of pair interactions of atoms 
in YBa2Cu307_5 crystal it can easily be shown that the 
abnormal temperature behavior of the elastic modulus 
C3333 is caused by change with temperature of force con- 
stants of the apical oxygen atoms. In used approach the 
modulus C3333 is determined by the interaction of apex 
oxygen atoms with copper and barium atoms but as ab- 
normal temperature behavior of the elastic moduli Cnn, 
C2222 is not observed, a sole opportunity is temperature 
renormalization of the 0(4) oxygen - copper force con- 
stants. The given conclusion is consistent with experi- 
mental results @T] (also references in it) [421 1431 14"4"| in 
which features in dynamics of O (4), 0(1) oxygen atoms 
in YBa2Cu307_5 were observed also. From the point 
of view of study of dynamics of atoms in YBa2Cti307_5 
the recently published article is of interest. In 
it is shown that in the temperature region 100 - 200 K 
incoherent lattice fluctuations essentially surpass ther- 
mal ones. Growth of incoherent lattice fluctuations u(T) 
were observed at temperature T* « 200.fr and the criti- 
cal temperature T* is attributed in to the formation 
of charge ordering. It is very interesting that the tem- 
perature interval in which incoherent lattice fluctuations 
exceeds thermal coincides with the interval in which hys- 
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teretic behavior of speed of ultrasound is observed. And 
the temperature of formation of charge ordering T* corre- 
sponds to the temperature of opening of a hysteresis loop. 
In this connection I shall pay attention to the experi- 
mental result 0Hj where hysteretic behavior of ultrasonic 
wave velocity in CuO crystals was observed. And recently 
in this compound charge-ordered domains were revealed 
|47| . Charge-ordered domains with necessity results in 
occurrence of correlated states of lattice degrees of free- 
dom, that may be the reason of formation in a crystal of 
strongly correlated groups of atoms. Recently, we offered 
the model of strongly correlated sublattice for an expla- 
nation of hysteretic behavior of ultrasonic wave velocity 
and thermal conductivity in YBa 2 Cu307_,5 crystals |5U] . 
Within the framework of the given model I have calcu- 
lated temperature dependence u(T) = yj a(T) (for more 
details see [SOI)- The calculated curve u(T) describes ex- 
perimental results for temperatures T < Tj*. However 
at higher temperatures on the experimental dependence 
u{T) essential growth of incoherent lattice fluctuations 
is observed. The given growth of fluctuations may be 
caused by a disordering of oxygen atoms in the Cu(l)- 
0(1) chains. Really, a quiet large anomaly of a linear 
expansion coefficient was observed in nearly optimally 
doped ((5=0.68) of YBa 2 Cu307_,5 untwined single crys- 
tals in the same temperature region [2Bj. But in fully 
oxygenated (5=7.0) crystals such anomaly were absent 
(see inset in fig. 3). I have revaluated the u{T) temper- 
ature dependence with taking into account a disordering 
of vacancies of atoms of oxygen (Fig. 3, curve 3). One 
can see that the calculated temperature dependence well 
enough describes experimental results. 

Thus, the analysis performed shows that ob- 
served experimentally incoherent lattice fluctuations in 
YBa2Cu307_,5 crystals may be caused by particular dy- 
namics of oxygen atoms and I can speculate that the 
hysteretic behavior of the C3333 elastic modulus may be 
attributed to the formation of charge ordering at tem- 
peratures T < T*. 

IV. Conclusion 

The joint analysis of temperature dependencies of 
sample resonant frequencies and ultrasonic wave velocity 
as in YBa2Cu307_,5 single crystal samples with various 
twinning structures and in polycrystal samples measured 
over a brode range of frequencies (100 kHz -=- 20 MHz), 
shows strong anisotropy of hysteretic behavior of elastic 
moduli. Namely, the C3333 and in much smaller degree 
C2323 moduli have hysteretic temperature dependence. 
While the moduli Cim, C2222, C1212 and C1313 do not 
show abnormal temperature dependence. The analysis 
of a tensor of elastic moduli of YBa2Cu307_,5 crystal 
on the basis of the microscopic model which is taking 
into account interaction of atoms of the nearest coordi- 



nation spheres has shown that the anomalies in behav- 
ior of an elastic modulus are caused by a temperature- 
dependent renormalization of force constants of apex oxy- 
gen atoms with Cu(l), Cu(2) copper atoms and appar- 
ently are attributed to the formation of charge ordering 
in YBa2Cu307_5 crystals. 
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